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MASSIVE YOUNG STELLAR OBJECT W42-MME: THE DISCOVERY OE AN INERARED JET USING 

VLT/NACO NEAR-INERARED IMAGES 

L. K. Dewangan\ Y. D. Mayya\ A. Luna\ and D. K. Ojha^ 

ABSTRACT 

We report on the discovery of an infrared jet from a deeply embedded infrared counterpart of 6.7- 
GHz methanol maser emission (MME) in W42 (i.e. W42-MME). We also investigate that W42-MME 
drives a parsec-scale H 2 outflow, with detection of bow shock feature at ^0.52 pc to the north. The 
inner ^0.4 pc part of the H 2 outflow has a position angle of ^18° and the position angle of ^40° is 
found farther away on either side of outflow from W42-MME. W42-MME is detected at wavelengths 
longer than 2.2 /im and is a massive young stellar object, with the estimated stellar mass of 19±4 Mq. 

We map the inner circumstellar environment of W42-MME using VLT/NACO adaptive optics and 
L' observations at resolutions ^072 and respectively. We discover a collimated jet in the inner 

4500 AU using the L' band, which contains prominent Bra line emission. The jet is located inside an 
envelope/cavity (extent ^10640 AU) that is tapered at both ends and is oriented along the north-south 
direction. Such observed morphology of outflow cavity around massive star is scarcely known and is 
very crucial for understanding the jet-outflow formation process in massive star formation. Along 
the flow axis, which is parallel to the previously known magnetic field, two blobs are found in both 
the NACO images at distances of ^11800 AU, located symmetrically from W42-MME. The observed 
W42-MME jet-outflow configuration can be used to constrain the jet launching and jet collimation 
models in massive star formation. 

Subject headings: ISM: jets and outflows - ISM: HII regions - ISM: magnetic fields - ISM: individual 
objects (W42-MME) - stars: formation 


1. INTRODUCTION 


The driving mechanism of jets and their role in massive 
star formation (MSE) is still subject to considerab le de¬ 
bate. Outflows are an nbiqnitons feature of MSE (Lada 




1985 Beuther et al. 


Zhang 2005). There have 


een increasing number of detecti on of molecular out 


flows from massive protostars (e. g. Beuther et al. 2002, 
2009 Wu et al.|2004[ Zhang|2QQ5 Arce et al.|2007p , that 


suggests the accretion process for their formation. How¬ 
ever, the number of known massive young stellar objects 
(M YSOs) associated with collimated jets is very limited 
(see [Guzman et ah 2010 for more details). Addition¬ 


ally, the information of magnetic field tow ards MYSO 


2004 

: Garrasco-Gonzalez et al. 

2010). The investigation 

of M 

Sb' is often handicapped 

by the fact that massive 


stars are statistically rare, situated at large distances 
(> 1 kpc), present in clustered environments, and suf¬ 
fer from high extinction. All these factors inhibit their 
observational study in the early phases. Key to overcome 
these difficulties is the high-angular resolution observa¬ 
tions towards candidate massive protostars at infrared 
and longer wavelengths, which will allow mapping of the 
inner circumstellar environment. It has been well estab¬ 
lished that the 6.7-GHz methanol maser emi ssion (MME) 


is associated with early ph ases of MSE (e.g. Walsh et al. 


1998 Urquhart et al. 2013) and hence traces the locations 
of MYSOs. Therefore, such sites are very promising to 
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carry out a detailed study to understand the involved 
physical processes of MSE. 

W42 is a Galactic giant Hll region (G25.38—0.18) 
harboring a spectroscopically classified 05 -06 star in 
the h eart of a near-infrared (NIR) cluster (Blum et al.| 


2000). Dist ances reported for W 42 range from 2.2 kpc 
obtained by |Blum et al.| (|2000|) assuming the 0-star is 
in the ZAMS, to 3.8 kpc obtained kinematically using 
radio re combination line H 86(a from the Hll region (59.1 
kms“^; Lester et al.|l985 ) , and GO line from the associ¬ 
ated cloud (58-69 kms“^; Anderson et al.||2009|. Using 
the K-b and polarimetric im age towards the NIR cluster 
in W42, Jones et al. (2004) discovered a bipolar reflec¬ 
tion nebula that is illuminated by an undetected source 
at ^10" south-west of the 0-star. Based on the align¬ 
ment of the long axis of the nebula with the magnetic 
field at the position angle of ^15°, they argued that the 
reflection nebula is tracing a bipolar outflow driven by 
an embedded young stellar object (YSO) that is hidden 
by a cloud with extinction of Ay > 55 mag. However, 
they did not find any K-band counterpart of the driv- 
ing source. We found that a 6.7 GHz MME detected by 


Szymezak et al. (2012) lies close to the predicted position 
of the driving source. We trace the reflection nebula as 
a parsec-scale bipolar outflow in the H 2 (1—^ S(l) 2.12 
/am continuum-subtracted image (see Eigure[H. The line 
of sight velocity of the MME is ^58.1 km s , very sim¬ 
ilar to the velocity of the Hll region. Gonsidering the 
fact that the 6.7-GHz maser is physically associated with 
W42, a distance of 3.8 kpc is also adopted for MME in 
this paper. The presence of a parsec-scale outflow around 
6.7-GHz MME makes this system ideal for exploring the 
early phases of MSE, including looking for a jet using 
high angular resolution NIR imaging. 
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In this paper, using multi-wavelength archival data, 
we investigate an infrared counterpart (IRc) of 6.7-GHz 
MME (hereafter W42-MME) and present observational 
evidence for W42-MME as a MYSO with a jet at the 
base of a parsec-scale bipolar H 2 outflow. The presence 
of jet is investigated using Very Large Telescope (VLT) 
NIR adaptive-optics images and the stellar counterpart 
is investigated using mid-infrared images. 


2. DATA AND ANALYSIS 
2.1. Adaptive-optics near infrared imaging data 

Adaptive-optics imaging data towards 6.7-GHz MME 
were retrieved from the ESO-Science Archive Eacilit}0 
(ESO proposal ID: 089.G-0455(A); PI: Joao Alves). The 
images were observed in Kg-band (Ac = 2.18/im, AA = 
0.35 jam) and L'-band (Ac = 3.80 jam, AA = 0.62 jam), us¬ 
ing 8.2m VLT with NAOS-GONIG A (NAGO) ada ptive- 
optics systenp] ( [Lenzen et al.||2003 Ronsset et al.||2003 ). 
Eive Ks frames and six L^ frames of 24 and 21 seconds, 
respectively, were used in this study. These data were 
reduced using the standard analysis procedure available 
in IRA E and STAR-L INK packages, as described in de¬ 
tail by Kumar ( 2013| ). The astrometry of NAGO images 


was cali brat ed using the GPS K-band point sources (see 
Section |2.2 for more details). The plate scales of final 
processed hi AGO and L' images were 07054/pixel and 
07027/pixel, respectively, resulting in a resolution of 072 
(^760 AU) and 07l (^380 AU), respectively. 


2.2. Other Archival Data 

We also utilized the multi-wavelength data from dif- 


tic R ing Survey (GRS: 


13< 


; Helfand et al. 

20061), Galac- 

CO(J = 1-0); 

Jackson et al. 

— —rTT—r^r—T- - 


(ATLASGAL: 870 jam; |Schuller et al.|j2009|, Herschel 
Infrared Galactic Pla ne Survey (HLGAL: 70-500 /im; 


Molinari et al. 2010), Galactic Legacy Infrared Mid- 
Plane Survey Extrao rdinaire (GLIMPSE: 3.6-8.0 jam 
Benjamin et al.||2003|), UKIRT Wide-held Infrared Snr 


vey for H2 (UWISHz: 2.12 /im; |Eroebrich et al. 2011), 
and UKIRT NIR Galactic P lane Survey (GPS: 1.25-2.2 


/im; Lawrence et al. 2007| )). To obtain a continuum- 


subtracted H 2 map (hereafter H 2 map), the GPS K-band 
image was registered and scaled to the UWISH2 narrow- 
band H 2 (1—0) S(l) 2.12 /im image (resolution ~1"). 
The H 2 map was smoothed with a Gaussian of sigma ^1 
pixel to enhance the faint featu res. The K-band p olari- 
metric data were obtained from Jones et al. (2004). 


3. RESULTS 


3.1. Multi-wavelength view around 6.7-GHz MME 

The distribution of molecular H 2 emission is shown in 
Eigure Ef. where the most prominent extended struc¬ 
ture is a continuous elongated emission feature centered 
on the position of the 6.7-GHz MME. On the northern 
side, the feature intercepts a bow shock at a projected 
distance of ^0.52 pc from the 6.7-GHz MME. The ex¬ 
planation of the origin of H 2 emission in W42 is difficult 
because an 05-06 star is located at ^0.22 pc distance 


^ http://archive.eso.org/eso/eso-archive_main.html 
^ http: / / www.eso.org/sci / facilities / par anal / instruments / naco. html 


from the MME in W42 H ll region. The presence of H 2 
emission can be interpreted by either ultraviolet (UV) 
fluorescence or shocks. Considering the observed mor¬ 
phology of the H 2 features, we suggest that the elongated 
H 2 emission is caused by shock, which traces a parsec- 
scale bipolar H 2 outflow. The presence of bow shock 
feature and knots embedded in diffuse emission along 
the elongated H 2 feature offers further evidence for the 
shock-excited outflow activity. The overall morphology 
of the H 2 feature resembles the bipolar outflows of Her- 
big Haro objects. The outflow nature of the observed 
H 2 feature is further supported by recent detection of 
[Ee H] emissi on coinciding wi th the southwest part (see 
Eigure 12 of Lee et al. 2014). The inner ^0.4 pc part 
of the H 2 outflow has a position angle of ^18 °, which 
coincides with the reflection nebula reported by [Jones et 
Z1( [2^ (see their Eigure 5). However, the II 2 feature 


continues much further than the reflection nebula on ei¬ 
ther side at a position angle of ^40°. Considering the 
variation in position angles along the outflow, the out¬ 
flow morphology appears a little twisted in “S” shape. 
The GRS ^^CO(J = 1—0) line data are unable to pro¬ 
vide any outflow signatures towards MME due to coarse 
beam (beam size ^45"). Additionally, there are no high- 
resolution CO observations available in the literature. 

In order to identify the IRc of 6.7-GHz MME, we 
present H 2 , 3.6 /im, and 5.8 jam images in Eigure [^. 
W42-MME is barely detected in the GPS K-band im¬ 
age (not shown here) with upper limit of ^12.815 mag. 
It is well detected in GLIMPSE 3.6-5.8 jam images but 
saturated in the 8.0 jam (see Eigures and [ip. The 
GLIMPSE-I catalog lists photometric magnitude only 
in one of the GLIMPSE bands. Therefore, aperture 
photometry was carrie d out for W42-MME using the 
GLIMPSE images (see Dewangan et al. 2015 for more 
details). The GLIMPSE photometric niagnitudes were 
measured to be 8.295T0.031 (3.6 /im), 6.178T0.012 (4.5 
/im), and 4.545T0.103 (5.8 jam) mag. The color excess 
crite ria further reveal this s ource to be a Class I YSO 
(see Dewangan et al.||2015[ for classification schemes). 
Figure also illustrates the distribution of ionized emis¬ 
sion using MAGPIS 20 cm map (beam size ^6"). It 
can be noticed that the MME-outflow system is im¬ 
mersed within the W42 Hll region. The GRS ^^CO(J = 
1—0) molecular, ATLASGAL dust continuum at 870 jam 
(beam size ^ 1972), and Herschel sub-millimetre emis¬ 
sions (beam size ^6-18") are also^peaked close to the 6.7- 


SQ^pe 

GHz MME position (see Figure]^. These data therefore 
suggest that W42-MME is located at the densest part of 
the molecular cloud in W42. 

Apart from the elongated bipolar feature, the H 2 map 
also shows H 2 emissions in the south-west and south di¬ 
rections, as highlighted by black arrows (see Figure [^). 
These particular H 2 features are nicely coincident with 
GLIMPSE infrared and 20 cm emissions. In general, 
GLIMPSE images (3.6-8.0 jam) have been used to trace a 
photodissociation region around the HII region through 
the presence of polycyclic aromatic hydrocarbon features 
at 3.3, 6.2, 7.7, and 8.6 jam. Therefore, it seems that 
these H 2 features (see black arrows in Figure f^) are 
formed due to UV fluorescence, tracing the wall of an 
ionized cavity-like structure (see Figure p!p). The ion¬ 
ized cavity-like structure seen in GLIML^E images is 
also shown in Figure Interestingly, the bipolar H 2 




































































Inner circumstellar environment of W42-MME MYSO 
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Figure 1. a) Continuum-subtracted H 2 (1—0) S(l) 2.12 /xm image 
illustrating the outflow, b) RGB color image (in linear scale) with an 
overlay of the 20 cm continuum emission (0.211 Jy/beam X (0.1, 0.2, 
0.3, 0.4, 0.55, 0.7, 0.85, 0.95)). In both the panels, the positions of a 
6.7-GHz MME (O) and an 05-06 star (x) are marked in the flgure. 
The 6.7-GHz MME appears at the center of a bipolar emission. The 
IRc coinciding with the 6.7-GHz MME (i.e. W42-MME) is clearly seen 
in the RGB image. The cyan arrows highlight noticeable features, in¬ 
cluding a bow shock, along a continuous elongated H 2 emission. Apart 
from the elongated H 2 emission, the black arrows also show H 2 emis¬ 
sions in the south-west and south directions. The long axis of the H 2 
nebulosity is aligned with the rn agnetic held {B) at the position angle 
of ~15° (e.g. [Jones et aL||2QQ4| ). The magnetic held direction is also 
shown by a tnick blue arrow. (Note: In the displayed H 2 map, the 
K-continuum is scaled so as to avoid negative values within the H ii 
region, which has resulted in the under-subtraction of stellar features.) 


feature in the south direction is further extended beyond 
the cavity-like structure (Figure]^. 

To estimate the stellar mass, we modelled the Spec¬ 
tral Energy Distribution (S EP) of W42-MME using an 
-line SEP modelling tool ( [Robitaille et al.|2QQ6 2007). 


on 



Figure 2. a) Overlay of the GRS ^^CO dotted contours on the 8.0 
p,m im ag e. The dot-dashed black box is shown as a zoomed-in view in 
Figure jd] b) Golor composite map using Herschel images (in logarith¬ 
mic scale). The ATLASGAL 870 /xm emissions are overlaid by orange 
contours with 10, 20, 30, 40, 55, 70, 85, and 95% of the peak value 
(i.e., 3.3 Jy/beam). The solid black box is shown as a zoomed-in view 
in Figure [^. In both the panels, the other marked symbols are similar 
to those shown in Figure^ The boundary of a cavity-like structure is 
shown by the 5.8 /xm cyan contour with a level of 555 MJy/sr in both 
the panels. 


Spitzer-GLIMFSE fluxes at 3.6, 4.5, and 5.8 /im bands 
were used as data points, whereas GPS K-band magni¬ 


tude and ATLA SGAL 870 /im flux (obtained from Gsen- 


geri et al. 


2014 ) were treated as upper limits (e.g. Dewan- 


gan et al. 2015). We provided the visual extinction in the 
range 0-70 mag and our adopted distance to W42-MME, 
as input parameters for SED modelling. The fitted SED 
models of W42-MME are shown in Figurej^. The stellar 
mass distribution for W42-MME is also plotted in Fig¬ 
ure It- The weighted mean values of the stellar mass 
and extinction are 19±4 Mq and 48T15 mag, respec¬ 
tively. The extinctio n value of W42-MM E is consistent 
with the estimates of Jones et al. (2004). The weighted 
mean value of age is found to be 8.1 x 10^ year. Our 
SED result suggests that W42-MME is a MYSO candi¬ 
date with a weighted mean luminosity of 4.5 x 10^ L0. 
We also fitted the SED assuming a distance of 2.2 kpc 
and estimated a weighted mean mass of 14±3 Mq. W42- 
MME is still a massive protostar candidate at a shorter 
distance of 2.2 kpc. The SED result of W42-MME is also 
in agreement with the generally accepted argument that 
the 6.7-G Hz methanol masers a ie solely associated with 
MSF (e.g. |Urquhart et al.||2013). 

Taken together, it appears that W42-MME is a driving 
source of a parsec-scale bipolar outflow and the axis of 
the outflow is parallel to the magnetic field (see Figure [^. 


3.2. Inner circumstellar environment of W42-MME 

In Figure we show the inner circumstellar envi¬ 
ronment of W42-MME using the VLT/NAGO adaptive- 
optics images at and L' bands. The image shows 
diffuse emissions near the position of 6.7-GHz MME, as 
well as, in north and south directions (see arrows in Fig¬ 
ure |^), which are coincident with the emissions seen in 
H 2 map. The L' image reveals a point source (i.e. W42- 
MME) associated with a small-scale feature (within a 
scale of ^4500 AU). The small-scale feature is notice¬ 
ably elongated in vertical direction and the full width 
half maximum of its point spread function is almost twice 
that of other stars seen in L' image. This comparison al- 
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Figure 3. a) The observed SEP (fille d symbols) for W42-MME, ob¬ 
tained with the |Robitaille et al.|(|2006|) fitting procedure. Only those 
models are selected that satisfy the condition (x^ - Xbest) data 
point < 3. The black solid line corresponds to the best fit; the gray 
solid curves represent all other models providing a good fit to the data. 
The dashed curve shows photospheric contribution, b) Distribution 
of the stellar mass for all plotted models in Figure We derive a 
weighted mean mass of 19±4 Mq. 


lows US to rule out the small-scale feature as a point-like 
source. W42-MME is also surrounded by an envelope¬ 
like emission, which is extended on a physical scale of 
^ 10640 AU, oriented a long the north-south direction. 


Zhang & Tan (2011) performed the radiative transfer 
calculations to understand the observed NIR and mid- 
infrared emissions around embedded MYSOs. They sug¬ 
gest that the outflow cavity is the most significant fea¬ 
tures in images up to 70 /im. They also find that the 
NIR emission is due mainly to scattered light from the 
outflow cavity wall. Using the prediction of model, we 
consider the envelope-like structure as an outflow cavity 
traced by the scattered light. The cavity has very narrow 
opening tip at both ends and is more extended (^4900 
AU in east-west) near the position of W42-MME (see 
Eigure|^). 

In both NACO images, the diffuse emissions are found 
parallel to the cavity in north and south directions, which 
are separated by a similar distance ^371 (^11800 AU) 
with respect to W42-MME. It provides the evidence for 
the presence of a symmetric configuration along flow axis. 
The physical association of small-scale feature, the cav¬ 
ity, and the outflow is shown in Eigurej^. In the north 
direction, the H 2 knot is coincident with the tip of the 
cavity along the flow axis. All these observed features 
suggest that they can be originated by the interaction of 
a jet with molecular gas in its surroundings. Therefore, 
the small-scale feature within the cavity is witness of a 
jet associated with MYSO W42-MME (see Eigurej^ for 
zoomed-in view). 

Considering together the observed emissions, NACO L' 
image traces the inner part of the outflow, which is highly 
collimated to a small-angle jet. A zoomed-in view of the 
cavity and the jet using a two color composite NACO 
image (red:L' and greemK^) is shown in Eigurej^. The 
IRc and the cavity have red colors suggesting absence of 
Ks emission from the parts. It also appears that the cav¬ 
ity resembles an onion-like structure. The jet is bright 
in L' image, which encompasses the Bra (4.05 /im) line 
emission as main contributor. Therefore, we suggest that 
it is an ionized jet associated with a driving source W42- 
MME. The diffuse emission is also coincident with 
jet-like feature located within the cavity, which could be 
due to the Bry (2.166 /im) line emission. In order to 
further investigate the ionized nature of the jet, we es¬ 
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Figure 4. VLT/NACO adaptive-optics images (in logarithmic gray 
scale) of W42-MME. The selected area of NACO images is shown by 
a box in Figure [^. a) Kg band. Diffuse emissions are highlighted by 
white arrows (al^ see Figure ^]d). b) Overlay of L' contours on the 
L' image (the outer most contour corresponds to 4cr, with successive 
contours increasing logarithmically). In both the panels, the position 
of the driving source (^2000 = 14® .54, ^2000 = —06°48^01^^86) 

is shown by a red circle. Distances of diffuse emissions seen in L' image 
are also highlighted with respect to W42-MME. 


timated the emission measure (EM) of jet feature using 
its photometric fluxes obtained in both the NACO im¬ 
ages. The photometric calculations of jet-feature were 
performed with the removal of its surrounding diffuse 
emission. In the calculation of EM, we assume that the 
whole jet emission is dominated by the hydrogen lines. 
Therefore the and L' fluxes of jet feature provide a 
measure of the Bry and Bra line emissions, respectively. 
Eollowing the recombination theory and the observed 
flux ratio of Bra to Bry, we computed the extinction 


at Br a and Bry emissions (see Appendix to Ho et al 


1990 and references therein). The intrinsic flux ratio of 


Bra to Bry was adopted f or Case B with Tp=lQ^ K and 


ne=I0^ cm Eollowing 


Indebetouw et al. 


(2005) ex¬ 


tinction law, the visual extinction (Ay) of jet emission 
is also estimated about 37 mag using the extinction at 
Bry emission, which is in agreement with our SED result. 
Eurthermore, the extinction corrected Bry flu x is used to 


measure the EM o f ^2.7 x 10^ cm ^ pc (see Tokunaga 
fc ThompsonjI97'9 for the equation), which can be taken 


as an indicative value. Our estimate of EM is consistent 
with the valu es obtained for candi dates radio jets toward 
MYSOs (e.g.|Guzman et al.||20I2|). It is also to be noted 
that the magnetic field is parallel to the jet-outflow axis 
(see Eigure[^). In the south direction, the blob emis¬ 
sions and the tip of the cavity are not perfectly aligned 
with the H 2 outflow axis. Eurthermore, the H 2 outflow 
has a little bend with respect to the jet axis. It further 
suggests the variation in position angles from small scale 
to large scale along the outflow axis. 


3.3. Summary and Discussion 

We studied the inner circumstellar environment of a 
deeply embedded IRc of 6.7-GHz MME in W42 (i.e. 
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Figure 5. a) Overlay of contours (in cyan color) on the continuum- 
subtracted H 2 image at 2.12 ^tm. The selected area of H 2 image is 
shown by a box in Figure [^. The L' contours are similar to those 
shown in Figure [^. The H 2 contours (in yellow) are superimposed in 
order to highlight faint H 2 features. The magnetic field (B) direction 
is also shown as Figurefm. b) Zoomed-in view of a two color composite 
NACO image (in logarithmic gray scale), which clearly illustrates the 
morphology of cavity and jet. 

W42-MME), as the driving source of a parsec-scale bipo¬ 
lar H 2 outflow. The parsec-scale H 2 outflow morphology 
appears a little twisted in “S” shape, with the position 
angle changing smoothly from ^18° near the IRc to ^40° 
at the ends. W42-MME is characterized as a MYSO 
with mass of 19±4 Mq. This is the first time we identi¬ 
fied W42-MME and characterized its stellar mass. The 
aim of this work was to understand the jet/outflow for¬ 
mation process in MSE. We discovered a collimated jet 
in the inner 4500 AU in W42-MME using the L' band 
to which Bra line emission contributes prominently, sug¬ 
gesting the jet is ionized in nature. The presence of jet 
was further confirmed by the detection of bow shock fea¬ 
ture and knots in diffuse emission in H 2 map along the 
flow axis (i.e. jet-related features). The jet is found 
inside a cavity (extent ^10640 AU) oriented along the 
north-south direction, which is tapered at both ends. 
Such a morphology appears to be responsible for colli¬ 
mating the jet/outflow. The diffuse nebular blobs seen 
in north and south directions are located at similar dis¬ 
tances (^11800 AU) from W42-MME, which can be con¬ 
sidered as a symmetric configuration along the flow axis. 
Previously published magnetic field direction is parallel 
to the jet-outflow flow axis. 

All these observed characteristics at a few thousands 
AU and a parsec-scale share some common features of 



ical mechanism for W42-MME jet-outflow system. 

The high angular resolution of NACO images around a 


MYSO W42-MME have enabled the detection of the in¬ 
ner jet-outflow configuration. Such a system is hardly 
known and is very crucial for understanding the jet- 
outflow formation process. We find that the morphology 
of outflow cavity in W42-MME is unique and rare, which 
is not seen at the inner environment of some well-studied 


lo w and massive stars (e.g. Valufia my et ^ Euller et 


al. 2001 Carrasco-Gonzalez et al.||2010p. Therefore, the 
interred spatial morphology of W42-MME can be used 
to constrain the jet launching and jet collimation models 
in MSE. 
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